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ABSTRACT 


Optimum  control  theory  is  applied  to  develop  a  guidance  law  for  homing 
missiles.  An  existing,  closed  form,  general  solution  of  the  "minimum  error 
regulator"  problem  is  applied  to  a  previously  solved  problem  which  uses  a  very 
simple  plant  model,  in  order  to  verify  its  applicability.  The  solution  method  is 
then  applied  to  a  system  that  includes  autopilot  lag,  and  in  this  case  the  optimum 
law  is  shown  to  differ  from  proportional  navigation. 
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Section  I.  INTRODUCTION 


The  technique  of  proportional  navigation  has  been  found  to  be  the  most 
satisfactory  method  of  guiding  homing  missiles  —  a  fact  established  by  engineers 
through  many  years  of  design  experience  rather  than  by  analytical  proof.  In 
proportional  navigation,  an  attempt  is  made  to  mechanize  the  following  equation: 

y  =  N’A  5  2  <  N’  5  , 

where  y  is  the  missile  angular  turning  rate  and  A  is  the  line-of-sight  angular 
rate  referred  to  inertial  space.  N’  is  denoted  as  the  "Navigation  Ratio,"  and 
again,  the  range  of  acceptable  values  has  been  developed  mainly  through 
experience.  It  is  interesting  to  note  that  recently  it  has  been  rigorously  deter¬ 
mined  by  means  of  modern  control  theory  that  proportional  navigation  is  indeed 
optimal  in  that  for  unconstrained  control  effort,  the  miss  distance  at  intercept 
is  minimized  in  the  mean-squared  sense.  The  correspondence  between 
proportional  navigation  and  optimum  control  has  been  demonstrated  by  Bryson, 
Ho,  and  Baron,1  Janus,2  and  Speyer,3  among  others. 

The  purpose  of  this  study  is  to  cast  the  simplified  homing  problem  solved 
by  Speyer  in  the  above  referenced  report  as  a  "minimum  error  regulator" 
problem.  This  problem  is  solved  in  closed  form  by  Ogata,4  who  used  the 
generalized  quadratic  index  as  a  performance  criterion.  It  will  be  shown  that  the 
Ogata  method  yields  identical  results  to  the  Speyer  variational  solution,  and  can 
thus  be  verified  to  be  applicable  to  the  problem.  The  problem  is  then  extended  to 
include  autopilot  lag,  and  a  new  solution  obtained,  from  which  several  interesting 
conclusions  can  be  drawn  regarding  the  form  of  the  control  law  and  the  nature  of 
the  navigation  ratio. 


'Y.  C.  Ho,  A.  E.  Bryson,  S.  Baron,  "Differential  Games  and  Optimal 
Pursuit-Evasion  Strategies,"  IEEE  Transactions  on  Automatic  Control.  AC-10, 
No.  4,  October  1965,  pp.  385-389. 

2J.  P.  Janus,  Homing  Guidance,  Aerospace  Corporation,  El  Segundo, 
California,  Report  No.  TOR-469(9990) -1,  December  1964. 

3J.  Speyer,  Optimal  Control  Theory  and  Biased  Proportional  Navigation, 
Raytheon  Corporation,  Bedford,  Massachusetts,  Memo  SAD-330,  November 
1967. 


4K.  Ogata,  State  Space  Analysis  of  Control  Systems,  Englewood  Cliffs, 
New  Jersey,  Prentice-Hall,  1964,  pp.  547-557. 
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Section  II.  THE  MINIMUM  ERROR  REGULATOR  PROBLEM 


Any  linear  dynamical  system  of  the  n  order  can  be  expressed  either  as 

th 

an  n  order  differential  equation  or  a  set  of  n  first  order  differential  equations. 
The  latter  is  known  as  the  state  formulation,"  and  is  used  herein,  since  it  lends 
itself  to  matrix-vector  notation  and  manipulation.  It  will  be  assumed  that  the 
system  differential  equation  is  given  by: 

x  =  Ax  +  Bu  ;  x(0)  =  C  , 


where 


x  =  n  dimensional  column  state  vector 
u  =  r  dimensional  control  vector 
A  =  n  x  n  matrix 
B  =  n  x  r  matrix 

and  where  the  following  index  is  to  be  minimized: 

T  T 

J(C,T)  =  x*  (T)Px(T)  +  /  x*  (t)Qx(t)dt  +  /  u*  (t)R(t)u(t)dt 

Terminal  State  Control  cost 

state  weighting 

weighting 

where  the  *  symbol  denotes  the  conjugate  transpose  of  the  vector,  or  simply  the 
transpose  for  real  vectors,  and  P,  Q(t),  R(t)  are  matrices  of  appropriate 
dimensions. 

Ogata  shows  that  the  optimum  controller  for  such  a  system  can  be 
obtained  by  solving  the  nonhomogeneous  matrix  Riccati  equation: 

ds  . 

-JjT  =  -sA  -  A*s  +  sBR-1  (0)  B*  s  =  Q(0)  . 

If  the  matrices  A  and  B  are  constant,  i.e. ,  if  the  system  is  stationary,  the 
above  matrix  equation  can  be  solved  in  closed  form  for  the  time  T: 
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B(T)  =  {[«ai(T)  +  022(T)P][0„(T)  +  4>a  (T) P I”1}  , 


FIGURE  1.  OPTIMAL  CONTROL  SYSTEM 
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The  solution  of  the  subject  problem  can  then  be  summarized  as  follows: 

7.)  From  system  state  equations,  and  the  given  performance  index,  the 
following  matrices  are  known: 

A,  B,  P,  Q,  R 


2)  The  matrix  M  can  be  formed  from  above 

-A  |  BR-!(0)B* 

M  = 

- 1 - 

Q(0)  |  A* 

MT 

3)  From  knowledge  of  M,  e  can  be  found.  There  are  several 

methods  for  computing  this,  but  the  use  of  the  Laplace  transform  is 
often  the  most  convenient.  This  method  involves  the  relationship 


«c[e 

where  s  is  the  Laplace  operator 
4)  Once  eMF  is  known, 

MT 


MT]  -  i.i  - 

and  I  the  unit 
all  of  the  <py 

*1I<T)  j 
*«<T>  ! 


M]-1 

matrix. 

are  obtained  from  the  relation 

012(T) 

022(T) 


From  this,  s(T)  can  be  computed: 


s(T)  =  { [ 02i (T)  +  022(T)P]  [0u(T)  +  0i2 (T) P]-1} 


5) 


Once  s(T)  is  known,  u  (t)  and  F(T  -  t)  are  immediately 
obtainable:  °‘3 


F(T  -  t)  =  R-’(t)B*s(T  -  t) 

Uopt(t)  =  F(T  "  t)X(t)  • 

The  optimum  system  thus  mechanized  will  minimize  the  given  quadratic  index 
in  the  solution  interval  0  ^  t  <  T.  It  should  be  noted  that  the  optimum  solu¬ 
tion  is  given  in  terms  of  time-to-go  (T  -  t)  rather  than  elapsed  time  t. 
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Section  III.  OPTIMUM  CONTROLLER  FOR  SIMPLIFIED  HOMING  SYSTEMS 


The  conventional  overall  homing  loop  can  be  depicted  as  shown  in 
Figure  2  and  can  be  readily  rearranged  with  target  acceleration  as  an  input  in 
the  manner  of  Figure -3-.  If  the  dynamic  lags  of  the  seeker,  autopilot,  and 
missile  are  neglected,  the  very  simple  model  of  Figure  4  is  obtained.  This 
figure  also  includes  an  exponential  decay  model  for  target  acceleration,  this 
being  the  model  used  by  Speyer. 5  With  reference  to  this  figure,  the  problem 
can  be  stated  as  follows:  Given  the  observable  states  y^,  y^,  andy’t>  the  con¬ 
trol  vector  is  determined  that  will  minimize  the  miss  distance  at  intercept 
only,  i.e. , 

yd(t)  !t  _  T  =  minimum  in  the  mean-squared  sense  subject  to  a  constraint 
on  available  control  effort. 


FIGURE  2.  TYPICAL  HOMING  BLOCK  DIAGRAM 

The  state  equations  can  be  written  directly  from  Figure  4: 

=  x2 

X2  =  gx3  -  gn^ 
x3  =  -2ux3 

5J.  Speyer,  op.  cit. 


1C 


FIGURE  4.  SIMPLIFIED  HOMING  SYSTEM 
or,  in  vector -matrix  form  x  =  Ax  +  Bu 
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0 

1 

0 

0 

x2 

- 

0 

0 

g 

*2 

+ 

-g 

*3 

0 

0 

~2u 

*3 

0 

The  general  form  of  the  minimum  error  regulator  law  is 

T 

J  =  x*(T)Px(T)  +  /  [x*  (t)Q(t)x(t)  +  u*  (t)R(t)u(t)  ]dt 

0 

and  must  be  adapted  to  the  specific  problem:  A  single  vector  xt  is  to  be 

minimized  at  intercept  only.  Also,  the  control  vector  has  a  single  component, 
and  the  control  constraint  matrix  reduces  to  a  scalar.  These  conditions  can 
be  achieved  by  defining  the  P,  Q,  R  matrices  as  follows: 


1 

0 

0 

0 

0 

0 

p  = 

0 

0 

0 

Q  = 

0 

0 

0 

0 

0 

1 

0 

_ 

0 

0 

0 

R  =  X  =  scalar; 

therefore 
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With  these  definitions,  the  index  J  becomes 


J  =  Xj2(T)  +  X  /  n  2dt  . 

0  1 

■j 

The  M  matrix  can  now  be  written,  since  all  of  its  components  are  known. 


0  -1 


0  -g 


0  2u 


0  gYx 


0  0  0  0 


IVx  X 

To  compute  e  ,  the  matrix  [si  -  M]  is  needed 


s  1 


0  s 


0  -g*A 


[si  -  M]  = 


0  0  s  -  2u  0 


0  0 


0  0 


s  +  2u 


This  matrix  must  now  be  inverted.  Since  it  is  almost  in  upper  triangular 
form,  the  simplest  method  of  inversion  is  to  append  a  unit  matrix  and  transform 
the  combined  matrix  to  upper  triangular  form  by  means  of  elementary  trans¬ 
formations.  This  step  yields 
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s 

1 

0 

0 

0 

0 

!  1 

1 

0 

0 

0 

0 

0 

0 

s 

g 

0 

-g 2 A 

0 

1  0 

1 

1 

0 

0 

0 

0 

0 

0 

s  -  2v 

0 

0 

0 

!  o 

1 

0 

1 

0 

0 

0 

0 

0 

0 

8 

0 

0 

j  0 

if 

0 

0 

1 

0 

0 

0 

0 

0 

0 

s 

0 

1 

1  0 

1 

0 

0 

Vs 

1 

0 

0 

0 

0 

0 

0 

s  +  2v 

I  0 

0 

0 

g/s3 

g/s 

1 

The  inverse  is  obtained  by  solving  the  above  matrix  as  six  equations  in 
six  unknowns,  using  each  of  the  six  column  vectors  to  the  right  of  the  dashed 
line  in  turn.  Since  the  coefficient  matrix  is  upper  triangular,  this  step  is 
rather  simple,  and  the  inverse  becomes: 


The  correctness  of  the  above  matrix  can  be  verified  by  observation  of  the 
product 


[si  -  M]_1[sl  -  M]  =  I  . 
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The  can  now  be  obtained  from  the  above  matrix  as  follows: 


(t)  =  is-1 


0  — 


0i2(t)  =  £"1y 


021  (t)  = 


022  (t)  =  £ 


-1 


_1_ 

S 

_1_ 

S2 

_S_ 


1_ 

L  g 

~  s2 

s2(s  - 

•J 

2u) 

1 

-g 

s 

s(s  - 

2u) 

0 

1 

s  -  2u 

i  1 

“-In, 

_jd 

s» 

0 

_2 

JL 

jsL 

0 

s3 

s2 

0 

0 

0 

0 

0  0 

0 

0  0 

0 

0  0 

1 

-k. 

S 

g 


s2(s  +  2u)  s(s  +  2v)  s  +  2u 


The  next  step  requires  computation  of  s(T  -  t) .  However,  since 
s(T  -  t)  =  s(t)  at  t  =  T  -  t,  we  may  compute  s(t)  first.  The  multiplications 
0i2P  and  022?  can  be  performed  before  taking  the  inverse  transform,  since  P  is 
a  number  matrix.. 
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The  a. . ,  k  in  the  above  matrices  are  symbolic  representations  for  the  inverse 
transform  of  the  actual  matrix  elements.  In  order  to  obtain  s(t) ,  the  matrix 
f^ii  +  ^i2p)  must  be  inverted.  Performance  of  this  step  on  the  symbolic 
representation  yields: 


-akb  baf  -  cae 
-a2k  aclc  -  a2f 
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Then  s  (t)  can  be  obtained  by  multiplication: 


akeh  -akbh  <baf  -  c'ae)h 


8(t)  ~  ak(ae*~~db)  akel  -akbl  (baf  *  cae)i 

y 

akej  -akbj  (baf  -  cae)j 

The  function  F(T  -  t)  can  now  be  readily  computed,  since 


F(T  -  t)  =  R“1B«s(T  -  t)  *— [0  -g  0]s (t) 


F(T  -  *>  =  Xak(ae~d~  '-akei8:  akblg  ‘<baf  '  cae)lgIt  =  T  -  t 


X(ae  -  bd) 


-e  b 


(bf  -  ce) 
k 


_ 1 1  =  T  -  t 


The  actual  vector  F(T  -  t)  can  be  obtained  by  evaluating  the  inverse  transform 
of  the  above  elements: 


•  ■  -i 


=  -t 


d  =  =  -&L 

a*  /  2X 


e  .  i  ,  i 


i  =  JB-1 


Introducing  above  expressions  into  the  F(T  -  t)  vector,  and  substituting  for  t: 


t  =  (T  -  t)  =  t 
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la  terms  of  the  original  system  parameters,  the  expression  for  u  ^  can  be 

written: 


opt 


+  N> 


-2vt 

e  g0  +  2ut  -  1 

J2 _ 


4u2t  2 

go 


n. 


The  above  results  are  identical  to  those  obtained  by  Speyer6  by  direct  use  of 
classical  variational  methods.  Thus,  the  two  techniques  are  equivalent. 


If  the  control  effort  constraint  is  removed,  i.  e. ,  X  is  made  zero,  the 
term  N'  becomes 


3g*t 


N»  =• 


g*t  2 
go 


3 


thus  reducing  to  a  conventional  navigation  ratio.  Also  since  the  angle  a  is 
approximately  equal  to  the  ratio 


A  = 


where  R  =  missile  to  target  range,  the  following  is  obtained: 

.  R*d  -  ya” 

A  - 

but  R  =  -V  where  V  is  the  closing  velocity,  and 
c  c 


Substitution  into  A  equation  yields: 
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or 


V  t  A  = 
c  go  ' 


This  can  be  substituted  Into  n.  equation  to  render 


N’V 

nr—rx  +  N’ 


-2ut 

e  80  +  2ut  -  1 
go 


4u2t  2 

go 


n. 


t 


which,  when  \  -  0,  is  recognized  to  be  the  conventional  proportional  navigation 
law  with  a  bias  correction  for  target  acceleration. 


The  results  to  now  are  not  new,  as  indicated  by  the  previously  noted 
references.  The  derivations  serve  as  an  independent  verification  of  prior 
work,  and  also  verify  the  adequacy  of  the  method  of  solution  used.  This 
mothod  will  be  applied  in  the  next  section  to  the  same  system,  augmented  by 
an  autopilot  lag. 
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Section  IV.  OPTIMUM  CONTROLLER  FOR  HOMING  SYSTEMS  WITH  AUTOPILOT  LAG 


Figure  5  depicts  the  same  simplified  homing  system  of  the  previous 
flection,  with  the  addition  of  a  single  time  constant  t  between  achieved  and 
commanded  accelerations.  The  state  equations  can  be  obtained  from  the 

figure  as: 


xi  =  x2 

*2  =  gX3  ~  gX4 


x3  =  -2uXg 


which  in  vector  matrix  form,  become 


• 

*1 

0  10  0 

Xi 

0 

• 

x2 

0  0  g  -g 

x2 

0 

= 

+ 

• 

*3 

0  0  -2u  0 

x3 

0 

• 

X4 

3 

i 

o 

o 

o 

x4 

Cl) 

-  — 

l—  -J 

—  _ 

mm  mm 

where 


A  1 

w  =  —  . 

T 


The  generalized  index  J 


T  . 

J  =  x*(T)Px(T)  +  /  [x*  (t) Q(t)x(t)  +  u*  (t)R(t)u(t)  ]dt 

0 


is  reduced  to  the  desired  form  by  defining 
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FIGURE  5.  AUGMENTED  HOMING  SYSTEM 


which  yields  the  same  form  as  before: 

T 

J  =  xt(T)2  +  \  /  n  2dt  . 

0  ° 


The  matrix  M  is  obtained  next: 


from  which  [si  -  M]  can  be  formed: 


-2u  0 


0  -co 
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This  matrix  is  almost  triangular  and  can  be  rather  easily  inverted,  by  the  method 
outlined  in  the  previous  section.  This  inverse  is  given  by 


19 


Correctness  of  this  matrix  can  be  verified  by  multiplying  it  by  [si  -  M]  and 
observing  that  the  unit  matrix  results.  The  0  are  obtained  as  before  from  the 
above  matrix,  since  J 

<MS)  012(S) 

[si  -  M]-1  = 

02l(s)  022  (®) 

The  next  step  in  computing  s(  t)  requires  formation  of  the  expressions 

[02l(t)  +  022  (t)  P  J 
[011  (t)  +  0i2(t)P] 


Since  P  is  a  number  matrix,  the  above  expression  can  be  formed  before  inverse 
transformation  of  the  0^: 


[02l(t)  +  022  (t)PJ 


0 

0 

0 

0 


0 

0 

0 

0 
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r*n(t)  +  0,a(t)P]  =£ 


S3(S-CU)(S  +  w)  +— •  -s2  s2(s  _2u)"  S2(S  “ 

g*  /  cuVx  \  ±  -g  _■  g  — 

-1  ”  s3  \(s  -  u)  (8  +  w)J  S  8  (s  -  2u)  8(8  -  0)) 


gV  _L 


/  o)2A 

\(S  -  w)  (S  + 


I  sz  \(s  -  w)  (s  +  w)/ 

Let  s(t)  of  the  above  matrix  be  denoted  symbolically: 


s  -  2u 

0  0 


1 

8-0) 


all 

a12 

a13 

a14 

a21 

a22 

a23 

a24 

0  0  a33  0 

a41  0  0  a44 

This  matrix  must  now  be  inverted,  since  s(t)  requires  (<Pa  +  012P)_1.  This 
inverse  is  shown  on  the  next  page,  where  the  A.^  and  B.  are  defined  as 
indicated: 


A  -  a22all  ~  a12a21 
1  &11 

A  =  a24all  ~  a21a14 
3  all 

as  =  -inis. 

5  all 

A  -  a44all  ~  a41a14 
7  all 

t>  _  A1A7  “  A5A3 

B‘ - X, 


A  =  a23aH  ~  a21a13 
2  ali 


a4  =  -■ 


Ac  —  - 


a41a13 

all 


A1A8  ~  A5A4 
A1 

A5A2  ~  A6A1 
Ala33 
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22 


The  correctness  of  this  inverse  is  extremely  difficult  to  check  by  the 
methods  previously  used.  However;  it  can  be  easily  checked  by  assuming  an 
arbitrary,  nonsingular  numerical  form  for  [<£n  +  ^P]  with  zeroes  in  the 
appropriate  positions.  The  inverse  can  then  be  computed  via  the  general  form 
of  page  22  and  then  multiplied  by  the  original  matrix.  If  the  unit  matrix  results, 
the  generalized  inverse  can  be  assumed  to  be  correct.  The  matrix 

12  11 

2  3  11 

0  0  10 

10  0  2 

is  of  the  required  form,  and  the  inverse,  as  computed  from  the  formula,  is 

4  2  1 

3  3  3 

_3_  2_  _1_ 

3  "3  "3 

0  3  0 

2_  1_  J_ 

3  "3  3 

which  is  readily  observed  to  be  correct. 

For  simplicity,  +  ^P)""1  will  be  denoted  by 

Fci  c2  c3  c71 


since  only  the  first  row  will  be  needed  for  s(t) .  The  latter  thus  becomes 


-2 

1 

0 

1 
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hC, 

he  2 

hC3 

hC4 

iCj 

1C2 

iC3 

iC4 

s(t)  » 

JC, 

jC2 

JC3 

JC4 

kCi 

kC2 

kC3 

kC4 

and  since 

F(T  -  t)  =  R^B+sfT  -  t)  =  *[0  0  0  u>]s(T  -  t)  , 

it  becomes 


F(T 


-  t)  -  — 


cokC  j 
X 


cokC2  cokC3  cokC  4 


_ It  = 


(T  -  t) 


The  optimum  controller  u  t  can  now  be  computed  from 


Uopt  =  -F(T  "  *> 


*1 


x2 


*4 


The  remaining  task  consists  of  obtaining  F(T  -  t)  in  terms  of  the  original 
system  parameters.  To  this  end,  the  gains  Ci  in  terms  of  the  a  are  first 
computed:  } 


£ _ a22a44 _ 

1  a44 (alla22  "  a12a2l)  "  a41  (a14a22  “  a12a24) 

*a12a44 

p  w  -  -  - 

a44  (a11a22  “  a12a2l)  “  a41  (a14a22  “  al2a24) 
c  _  a4</a33  ( a12a23  ~  a22aia) _ 

3  a44  (a11a22  "  a12a2l)  “  a41  C a14a22  “  a12a24 ) 

q  _ _  a12a24  ~  a22a14 _ 

4  a44(atia22  -  a12a21)  -  a4i(ai4a22  -  aj2a24) 
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25 
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OJt 


The  common  e  terms  cancel.  Additionally,  if  the  expression  is 
multiplied  and  divided  by  gt  2,  the  solution  form  becomes: 


where 


Ft  =  -N’ 
S° 


1  1  J_  1 2utgo6 


2ut  2ut 

S°  -  e  g°  +  1 


gt  2  gt  4v2 
6  go  B  go 


cot  cot 


2ut 


1  fcot  e  g0  -  e  ®°  +  1 

_  —I  go _ 

"  CO2  \  cot 

e  «°t  2 

go 


g°t  2 

go 


-6co3g2t  2  +  6co2g*t  2  1 1  -  e 


-cot 


go 


-N'  = 


go 


go 


(I) 


and  (I)  is  as  previously  defined. 

The  optimum  controller  u  is  computed  from 


=  n 

c 


Therefore,  in  terms  of  the  original  system  parameters,  the  control  vector 
becomes: 


u  =  -Ft 
opt  go! 


*1 

*2 

*3 

X4 
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but  from  the  previous  section,  it  was  shown  that: 


9 


V  t  a 
c  go 


9 


and  therefore  n  becomes 
c 


The  first  three  terms  are  identical  to  those  previously  obtained  for  the 
non-augmented  case  and  represent  biased  proportional  navigation.  The  last 
term  is  an  additional  bias  term  due  to  the  autopilot  lag,  and  it  is  interesting  to 
note  that  it  is  identical  in  form  to  the  target  bias  term.  It  should  be  noted  that 
when  the  lag  is  removed  (i.e. ,  r  —  0  or  w  —  00 ) ,  the  last  term  drops  out, 
since 
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The  nature  of  N',  the  navigation  ratio  bears  further  study.  This  quantity 
was  defined  on  page  27.  When  divided  by  2cu3,  it  can  be  rewritten 


Nf  = 


3g h  2/  -cot  \ 

3g*t  3 - ZjJ2.i1  -'e  «° 

_  go  a)  \ _ / _ 

3g2t  6g*t  e  3g2t  2  _  9  0  9 

3\  +  g*t  «  ...  *°.  -  —g - S2_+i£L-i£?£ - 


If  the  autopilot  lag  is  eliminated  by  letting  co  —  80  (r  —  0) ,  the  above 
reduces  to: 


N'  =■ 


3g2t  3 
*  go 

3\  +  g2t  3 
go 


which  is  identical  to  the  results  previously  obtained.  Thus,  the  entire  optimum 
controller  reduces  to  the  previously  derived  quantities  when  the  lag  term  is 
eliminated.  Where  to  is  finite,  however,  it  seems  to  play  a  significant  part  in 
the  value  of  N\  In  order  to  observe  the  behavior  of  N'  as  a  function  of  time 
and  co,  let  \  =  0,  i.e. ,  remove  the  constraint  on  control  effort.  If  this  is  done, 
the  g2  terms  drop  out,  and  N’  becomes 


N'  = 


3t  2  / 

.  3 

-  g0  1  -  e 

go 

O)  ' 

-cut _ 

6t 

e  g  3t 

go  go 


Igo  +  JL 

co  2uc 


-2cot 

3e  g0 


For  very  large  values  of  t  ,  the  above  converges  to 


3t  3 

N'  = — i2_ 

t  3 

go 


For  typical  values  of  t  ,  the  above  equation  was  programmed  with  several 

go 

values  of  co  used  as  a  parameter.  The  resulting  plots  of  N’  versus  t  are  shown 

go 

in  Figure  6,  from  which  one  concludes  that  for  significant  autopilot  time  con¬ 
stants,  N'  deviates  considerably  from  the  limit  value  of  3,  and  is  also 
nonstationary. 
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NAVIGATION  RATIO  (M*) 


FIGURE  6.  PLOTS  OF  NAVIGATION  RATIO  VERSUS  TIME  TO  GO 
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Section  V.  CONCLUSIONS 


It  has  been  shown  herein  that  the  Ogata  Sjolution  of  the  minimum  error 
regulator  problem  can  be  used  to  solve  the  homing  problem.  This  was  done  by 
independently  solving  a  previously  solved  problem  and  obtaining  identical  results. 
The  subject  method  of  solution  is  (at  least  in  principle)  straightforward.  For 
high  order  systems,  the  computational  aspects  become  menacing,  but  could  be 
handled  by  use  of  the  digital  computer. 

The  effect  of  autopilot  lag  was  also  studied  using  the  subject  method,  and 
several  significant  conclusions  can  be  reached: 

1)  The  navigation  gain  is  a  nonstationary  function  that  varies 
significantly  with  the  autopilot  lag. 

2)  The  optimum  controller  requires  an  additional  bias  term  appended 
to  the  conventional  biased  proportional  navigation  vector.  This 
term  is  Identical  in  form  to  that  resulting  from  target  acceleration. 

3)  The  added  term  requires  that  the  actual  missile  acceleration  be 
measured,  since  this  is  a  required  state.  However,  in  cases  where 
this  state  is  not  explicitly  measured,  the  commanded  acceleration 
could  be  passed  through  a  synthetic  time  constant.  The  latter  would 
be  tailored  to  match  autopilot  response,  and  thus  an  approximate 
measure  of  could  be  obtained. 


31 


UNCLASSIFIED 


DOCUMINT  CONTROL  DATA  RAD 

ft— Mffiy  •/  «/  a^iliact  ai^  /mtoxfni  —  tofl—  amt*  b#  —  f— ¥  rMw  t—  ovrmll  npott  !•  cfa— //i»Q 


— .  RIPORT  BKCUniTV  CLASSIFICATION 

Unclassified 


i-1  r*j  •}  •  rf* 


» 


Redstone  Arsenal .  Alabama  35809 


I.  RIBOAT  Tt  TLK 


OPTIMAL  CONTROLLERS  FOR  HOMING  MISSILES 


*■  OKtCNtPTIVC  NO' 

None 


AU  THONIS)  (Fittl  I 

G.  Willems 


•I  md  iMhMJn  daft) 


•  f  IIMW) 


?•.  TOTAL  NO.  OF  I 


1 7b.  NO.  OF  RKFI 


CONTRACT  OR  BRANT  NO. 


B*.  ORIGINATOR**  RIRORT  NUbflBCRtJ) 


a  onoj.ct  NO^DA)  1X263302D212 
AMC  Management  Structure  Code  No. 

**  5231.12.171 


RE-TR-68-15 


10.  D1GTR1  RUT  (ON  ITATCMCNT 

This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign  govern¬ 
ments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  this  Command, 

ATTN:  AMSMI-RE. 


12.  SPONSORING  MILITARY  ACTIVITY 


Same  as  No.  1 


Optimum  control  theory  is  applied  to  develop  a  guidance  law  for  homing  missiles.  An 
existing,  closed  form,  general  solution  of  the  "minimum  error  regulator"  problem  is  applied 
to  a  previously  solved  problem  which  uses  a  very  simple  plant  model,  in  order  to  verify  its 
applicability.  The  solution  method  is  then  applied  to  a  system  that  includes  autopilot  lag,  and 
in  this  case  the  optimum  law  is  shown  to  differ  from  proportional  navigation. 


-1473 


RCPLACBS  DO  FORM  14TB.  1  JAN  #4,  WHICH  IS 
OMOLCTI  FOR  ARMY  USB. 


UNCLASSIFIED 

lwcurity  sin  cation 


